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Abstract 

I discuss the problem of phase separation in cuprates from the point of view of 
the Landau theory of Fermi hquids. I calculate the rate of growth of unstable 
regions for the hydrodymanics and collisionless limit and, in presence of long 
range Coulomb interactions, the size of these regions. These are analytic 
results valid for any strength of the Landau parameters. 
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There is experimental evidence that in the cuprates there exists a phase separation 
between hole rich and hole poor regions Also numerical studies in the t — J model 
(which is believed to describe the low temperature physics of these materials) show clearly 
signs of phase separation 0. It is well established that at zero doping and temperatures 
higher than the critical temperature the cuprates can be well described as a Heisenberg 
antiferromagnet. This has been confirmed theoretically from the studies of non-linear sigma 
models (which reproduces the long wavelength behavior of the Heisenberg model) and the 
scaling laws of many measurable quantities . At finite doping holes start to move in the 
antiferromagnetic background and interact by exchanging spin waves. At the level of the 
t — J model the physics at long wavelengths is described by the Shraiman-Siggia model 0. 
In the last few years Emery and Kivelson [|l| have also studied the problems of holes in an 
antiferromagnetic background and show that the phenomena of phase separation is quite 
natural. 

In this paper I study the physics of phase separation using the tools of the Landau 
theory of Fermi liquids. This work is inspired by the problem of growth of instabilities in 
three dimensional Fermi liquids studied by Pethick and Ravenhall [0. Thus, I will make the 
assumption of a Fermi liquid behavior for the holes. This is actually not a strong assumption 
because I am going to use a part of the Landau theory which depends only on conservation 
laws for the distortions of the Fermi surface. Recently it was shown that this part of the 
Landau theory can also be applied in Luttinger liquids in integrable systems [^. In any 
event, analytic studies of the ferromagnetic region of the Hubbard model at strong coupling 
(close to the antiferromagnetic instability) have shown that the holes form a Fermi liquid 

i- 
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I will consider a two dimensional system of holes (spinless) with density n = ^ where pp 
is the Fermi momentum (n is the doping density). Since I am interested in the low doping 
regime I will assume that the Fermi surface is spherical and I will not take into account 
lattice effects. The density of states at the Fermi level is given by A^(0) = ^ where m* is 
the effective mass of the quasi-holes and m is the bare mass. 
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I will consider angle depend distortions of the Fermi surface (in this case a line) and 
parametrize the distortions by an angle 6, 

oo 

Sn{9) = uq + 2^^Un cos{n6) , (1) 

n=l 

where I expanded the deviations in Fourier series with coefficients m„. Also, it is assumed 
that the quasi-holes interact through the Landau parameters f^^f which in principle can 
be obtained from the experiment or can be calculated from some theoretical model (for 
instance, the t — J model). It is also useful to consider the dimensionless Landau parameters 
= iV(0)/p,p.. 

Since the interaction between the particles at the Fermi surface only depend on the 
angle between the momenta of the quasi-holes we can also expand the Landau parameters 
in Fourier series, 

oo 

Fj^^j, = Fo + 2 5] F„ cos(n(e - 9')) (2) 

n=l 

where Fn are the Fourier components. 

At large densities a fermionic liquid is stable against phase separation and density fluc- 
tuations which results in a positive compressibility. As we decrease the density it can enter 
a metastable region where it separates into a liquid-gas phase but it is still stable against 
density fluctuations of long wavelengths. In this phase the compressibility is positive. As the 
compressibility approaches zero (the so-called spinodal line) the liquid becomes unstable to 
long wavelength density fluctuations as well. In the context of the Fermi liquid theory this 
result can be understood from the expression for the sound velocity, c^, in a Fermi liquid 
which is given by 0, 

2 _ n(l + Fq) _ m* I + Fq 2 
~ mN{0) ~ m 2 ^'^^ 

where fi? = ^ is the Fermi velocity. Observe that if Fq < — 1 then c| < and an unstable 
collective mode grows in the system and there is a phase separation. For the moment I am 
considering as a starting point a homogeneous system and local interactions. It means that 
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the instability grows in the system as a whole. If we include the Coulomb interaction, as 
I will show later, then the instability has a finite size. For the moment we consider local 
interactions since the inclusion of Coulomb interactions is trivial. 

It is also well established that in the cuprates the effective mass of the holes can be 
many times the bare mass. In the context of the Fermi liquid theory this is due to the a 
non-isotropic interaction on the Fermi surface. Indeed, the relation between the effective 
mass and the bare mass is given by 

m* ^ f ^''^^i.P') P ' 
m ^ ^'^ de{p') 

where 'nP{p) = 0(/i — e(p)) is the zero temperature distribution for the holes and /i is the 
chemical potential. Substituting (2) and using the orthogonality relation between the Fourier 
components one finds, 

— = 1 + Fi. (4) 

m 

Therefore, the change in the effective mass depends only on the second Fourier component. 
From now on we will keep only Fq (related to the instability) and Fi (related to the increase 
in the mass) since these are the main ingredients for the physics of the system. 

Observe that for Fq < —1 the hydrodynamical rate of growth of the instability in the 
system (which determines the time scale for the growth, F"^) is obtained straight from (3), 

( [I + F,)\l + F,\ \"' 

r(g) = I ^ 1 qvF- (5) 

This result is valid when the collisions dominate the physics, that is, the mean free path is 
short. However, since the density of holes is small the mean free path can be quite large 
and the collisionless limit is more interesting. 

The collisionless limit can be obtained from the Landau equation for sound waves in the 
absence of a colhsion integral (which can also be included 0]). It reads 0, 

+ V, . V6n, - -vlu + Y. f,.^Sn)i = 0, (6) 



where U is an external field and in terms of the deviations (1) we have Sn^ = (~ ^ae(p)^ ) ■ 
Assuming an harmonic evolution in time, that is, 6n{9) ~ e**^^''~^*\ the Landau equation 
has the form, 

{qvF cos(^) - u) 5n{d) + qvp cos{0) i N{0)U + ^^(^ " d')5n{d') J = 0, (7) 
which can be simplified by our assumption of keeping only the two first Fourier components, 
{qvp cos{e) - u) 6n{e) + qvp cos(^) {N{0)U + FqUq + 2FiUi cos(0)) = 0. (8) 
Integrating (8) over 9 one finds, 

which is essentially the conservation of the number of particles. 
Solving (8) for 6n{9) one gets, 

^ qvp cos(g) ^ ^ cos(^)) , (10) 

uj — qvpcos{0) 

and an integration over 9 leads to, 

^^0 = Xrpa{s) {N{0)U + Fouo + 2sFiUi) (11) 

where s = and 

qVF 

, , /■2'^ (19 cos(9) s ^ 

XRPAis) = / 4- = -^== - 1 12 

Jo 2ns — cos{9) ^J — 1 
is the RPA density-density correlation function [Q. 

Since we are interested in the growth of collective modes we calculate the density- density 
correlation function for the problem, 

/ X _ ^0 Xrpa{s) , . 

Nm l-xnPA{s){Fo + s'A,) ^ > 

where. 



As it is well known the collective modes are given by the poles of the correlation function 
(13), that is, the points where = 0. Unlike the three dimensional case 0, which leads 
to a transcendental equation that can be solved only numerically, in two dimensions we can 
solve this problem analytically for any strength of the interaction. The equation has two 
roots which read, 



2 _ -1 + 2v4i - 2Fo + 2AiFo + AA^ + 4Fo + AF^ + AA^Fq 
'± " 2A,{2-A,) ■ ^^^^ 

It is easy to conclude that > and therefore is related to a true collective mode in the 
system and < which means that the rate of growth of the instability in given by, 

T{q) = s.qvp. (16) 

The inclusion of Coulomb interaction, V{r) = ^ (where e is the dielectric constant), is 
rather trivial since the interaction is isotropic. The Fourier transform of this interactions in 
two dimensions is given by V{q) = and following the prescription of the Landau theory 
0] the only change in our calculation would be to make the substitution, Fq Fq + 
Observe that this modification would lead to a momentum depend relation for the build up 
of the instability {Fq + ^^'^^^W < — 1) namely, 

e|i + -^ol 

This relation means that the instability can be probed only for wavelengths smaller than 
2^^^^^°i N which established the minimum size for the unstable regions in the plane, that is, 

e^m(l+Fi) ft'! 

the size of the role rich regions. Of course, in the absence of Coulomb forces (e = 0) the 
instability region has the size of the system. 

In conclusion, in this paper I use the Landau theory of Fermi liquids in order to calculate 
analytically the rate of growth of phase separation in a two dimensional systems of holes in 
cuprates for any strength of the Landau parameters. I also estimate the size of the unstable 
regions due to the Coulomb repulsion between the holes. 
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